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STEADY  STATES  AND  DYNAMICS  OF  2-D  NEMATIC  POLYMERS 
DRIVEN  BY  AN  IMPOSED  WEAK  SHEAR* 

HONG  ZHOUt  AND  HONGYUN  WANG* 

Abstract.  We  study  the  2-D  Smoluchowski  equation  governing  the  evolution  of  orientational 
distribution  of  rodlike  molecules  under  an  imposed  weak  shear.  We  first  recover  the  well-known 
isotropic-to-nematic  phase  transition  result  [G.  Marrucci  and  P.  L.  Maffettone,  Description  of  the 
liquid- crystalline  phase  of  rodlike  polymers  at  high  shear  rates,  Macromolecules,  22,  4446-4451,  1989]: 
in  the  absence  of  flow  the  isotropic-nematic  phase  transition  occurs  at  U  =  2  where  U  is  the  normalized 
polymer  concentration,  representing  the  intensity  of  the  Maier-Saupe  interaction  potential.  Then  we 
show  that  in  the  presence  of  an  imposed  weak  shear  there  is  a  threshold  (C/q  ~  2.41144646)  for  U: 
When  U  <Uo,  steady  state  solution  exists;  otherwise  there  is  no  steady  state.  Furthermore,  we  carry 
out  multi-scale  asymptotic  anlaysis  to  study  the  slow  time  evolution  driven  by  the  weak  shear.  It  is 
revealed  that,  to  the  leading  order,  the  order  parameter  of  the  orientational  distribution  is  invariant 
with  respect  to  time  whereas  the  angular  velocity  of  the  director  is  position-dependent.  When  U  <  Uq-, 
the  director  of  the  orientational  distribution  converges  to  a  stable  steady  state  position;  when  U  >Uq, 
the  angular  velocity  of  the  director  is  always  positive  and  the  orientational  distribution  will  not  reach 
a  steady  state.  Finally,  the  effect  of  weak  shear  on  the  phase  diagram  is  investigated.  It  is  found  that 
the  phase  relation  under  weak  shear  can  be  obtained  from  the  pure  nematic  phase  relation  through 
a  simple  algebraic  transformation. 

Key  words.  Nematic  polymers,  Kinetic  theory,  Maier-Saupe  potential.  Shear  flow,  Perturbation 
analysis 
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1.  Introduction 

The  rheological  properties  of  nematic  liquid  crystals  are  strongly  affected  by  the 
dynamic  behavior  of  the  molecular  alignment  [1,  20].  The  standard  Doi-Hess  kinetic 
theory  treats  the  nematogenic  molecules  as  rigid  rods,  and  describes  the  ensemble  with 
an  orientational  distribution  function.  The  dynamics  of  the  orientational  distribution 
function  is  governed  by  the  nonlinear  Fokker-Planck  (Smoluchowski)  equation.  The 
steric  mutual  interaction  of  rodlike  molecules  is  modeled  by  the  Maier-Saupe  interac¬ 
tion  potential  [26,  27,  4].  The  complex  behavior  of  3-D  nematic  polymers  subjected 
to  an  imposed  shear  flow  has  been  investigated  numerically  by  solving  the  Smolu¬ 
chowski  equation  using  spherical  harmonic  expansions  [13,  14,  15,  16,  17,  18,  19,  22]. 
However,  theoretical  analysis  on  this  subject  is  still  at  its  early  stage.  Recently  the 
mathematical  society  started  to  investigate  the  Smoluchowski  equation  theoretically 
[5,  6,  7,  8,  10,  11,  12,  24,  25,  29,  30,  31,  32].  The  purpose  of  this  work  is  to  extend 
the  analytical  studies  to  Smoluchowski  equation  under  an  imposed  weak  shear. 

The  paper  is  organized  as  follows.  Section  2  is  a  brief  review  of  the  results  for 
pure  nematic  polymers.  In  section  3  we  study  steady  states  of  nematic  polymers  in 
the  presence  of  an  imposed  weak  shear.  Section  4  covers  the  slow  time  evolution  of 
orientational  distribution  using  the  approach  of  multi-scale  asymptotic  analysis.  In 
section  5  we  investigate  the  behavior  of  the  Smoluchowski  equation  near  the  isotropic- 
nematic  phase  transition  (U  =  2).  The  main  results  of  the  paper  are  summarized  in 
section  6. 
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2-D  NEMATIC  POLYMERS  UNDER  AN  IMPOSED  SHEAR 


2.  Pure  nematic  polymers 

For  polymer  rods  in  the  two-dimensional  (2-D)  space,  the  orientation  of  a  polymer 
rod  is  described  by  an  angle  6  and  the  orientation  direction  is  mathematically  repre¬ 
sented  by  the  unit  vector:  u=  (cosd,sin0).  The  Maier-Saupe  interaction  potential  is 
given  by  [26]: 


Vms{0)  =  -ksTU  [{cos2e)cos2e+{sm20)sm2e] ,  (2.1) 


where  U  is  the  normalized  polymer  concentration,  representing  the  strength  of  the 
Maier-Saupe  interaction  potential,  and  ksT  is  the  Boltzmann  constant  times  the 
absolute  temperature.  For  mathematical  convenience,  in  this  study  we  set  ksT  =  1  (or 
equivalently,  all  potentials  are  normalized  by  ksT).  In  the  Maier-Saupe  interaction 
potential  (2.1),  (cos20)  denotes  the  average  of  cos20  with  respect  to  the  current 
orientational  distribution  of  polymer  rods  f{0,t).  The  time  evolution  of  f{0,t)  is 
governed  by  the  Smoluchowski  equation  [9] : 


df{0,t)  d_ 
dt  ''80 


fdV{0) 

80 


80  ) 


(2.2) 


where  Dr  is  the  rotational  diffusion  coefficient  of  polymer  rods  and  V (0)  the  sum  of  all 
potentials  affecting  polymer  rods.  If  the  Maier-Saupe  potential  is  the  only  potential 
affecting  polymer  rods,  then  at  equilibrium  the  orientation  distributation  of  polymer 
rods  is  given  by  the  Boltzmann  distribution: 


/(6»)~exp(-VMs(6'))- 


(2.3) 


Notice  that  when  the  Maier-Saupe  potential  is  the  only  potential  affecting  polymer 
rods  and  when  only  the  equilibrium  distribution  is  concerned  we  can  select  a  suitable 
coordinate  system  such  that  (sin26l)  =  0.  So  in  this  special  coordinate  system  the 
Maier-Saupe  potential  (2.1)  is  simplified  to 

Vms{0)  =  -U{cos20)cos20.  (2.4) 


In  this  special  coordinate  system,  the  Maier-Saupe  potential  is  completely  specified 
by  a  parameter  r  in  the  form  of  Vm3{0)  =  — rcos20  where  r  is  defined  as 

r=U{cos20).  (2.5) 

In  fact,  s=  (cos 20)  is  the  Flory  order  parameter,  where  s  =  0  corresponds  to  isotropic 
phase  and  s  =  l  describes  the  fully  aligned  phase.  The  equilibrium  distribution  (2.3) 
becomes 


m=  .  (2.6) 

/p  exp(rcos20)(i0 

The  probability  density  f{0)  given  above  depends  on  r.  Thus,  although  (2.5)  is  the 
definition  of  r,  it  is  also  an  equation  for  r  in  which  r  appears  on  the  left  side  as  well 
as  on  the  right  side  through  the  probability  density  f{0).  Rewriting  equation  (2.5)  in 
the  form  of  the  order  parameter,  we  have 


r 


U' 


(cos  20) 


(2.7) 
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Applying  integration  by  parts  in  the  numerator  of  (cos  20),  one  has 

1 


^27r 


^27r 


cos20exp(rcos20)d0 = 


exp(rcos20)dsin20 


(.2ir 


=  r  sin29exp{rcos20)d9. 
Jo 


Substituting  into  equation  (2.7)  and  dividing  by  r  yields 

(sin^20)  =  ^. 


(2.8) 


As  an  equation  for  r,  (2.8)  is  equivalent  to  (2.5)  when  r^O.  Equation  (2.8)  has  the 
advantage  that  r  appears  only  on  the  left  side  of  equation.  It  is  straightward  to  verify 
that  (sin^20)j.=o  =  5-  In  [7)  32],  it  was  proved  that 


dr 


(sin2  20)<O  for  r>0. 


(2.9) 


Therefore,  equation  (2.8)  has  a  nematic  solution  (r  >  0)  if  and  only  if  1/  >  2  [7,  25].  For 
t/>2,  there  is  a  unique  function  r{U)  such  that  equation  (2.8)  is  satisfied.  Further¬ 
more,  r(U)  is  an  increasing  function  of  U.  This  one-to-one  correspondence  between 
U  and  r  also  allows  us  to  view  U  as  a,  function  of  r: 


U{r)  = 


or  equivalently  U{r)  = 


1 


(2.10) 


(cos2^)  ^s^n^26')' 

The  right-hand  side  of  (2.10)  depends  on  r  only.  So  U{r)  can  be  calculated  in  a 


1.5 


2.5 


3.5 


Fig.  2.1.  Phase  diagrams  for  the  2D  pure  nematic  polymers  (no  shear  flow).  The  left  panel 
shows  r  as  a  function  of  the  normalized  polymer  concentration  U  where  r  is  defined  as  r  =  U (cos26) . 
Mathematically,  it  is  more  convenient  to  treat  r  as  the  independent  variable  and  treat  U  as  a  function 
of  r.  The  right  panel  shows  the  order  parameter  s  as  a  function  of  U  where  the  order  parameter  is 
defined  as  s  =  (cos20).  The  isotropic-nematic  phase  transition  occurs  at  U  =  2. 


straightforward  way  without  solving  non-linear  equations.  This  illustrates  that  it  is 
mathematically  more  convenient  if  we  use  r  as  the  independent  variable  although  r 
does  not  have  a  clear  physical  meaning  (it  is  the  product  of  the  normalized  poly¬ 
mer  concentration  and  the  order  parameter).  Figure  2.1  A  shows  r  as  a  function  of 
U  whereas  Figure  2. IB  depicts  the  order  parameter  s=(cos20)  as  a  function  of  U. 
Clearly,  the  isotropic-nematic  phase  transition  occurs  at  U  =  2. 
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3.  Nematic  polymers  under  an  imposed  weak  shear 

In  the  presence  of  an  imposed  2-D  shear  flow  v  =  j(0,—x),  where  7  is  the  shear 
rate,  the  non-dimensional  Peclet  number  is  defined  as  £  =  ^  where  Vr  is  the  rota¬ 
tional  diffusion  coefficient.  Here  we  use  £  to  denote  the  Peclet  number  because  we 
study  the  case  of  small  Peclet  number  in  this  paper.  Under  the  imposed  shear,  the 
orientation  distribution  f{0,t)  evolves  according  to  the  Smoluchowski  equation  [9]: 

+  +  •  (3.1) 

Recall  that  ksT  has  been  set  to  one.  The  rotational  diffusion  coefficient  Dr  can  be 
absorbed  into  a  non-dimensional  time  t  =  tDr.  For  mathematical  simplicity,  we  set 
Dr  =  l  (or  equivalently  we  use  the  non-dimensional  time).  The  Maier-Saupe  inter¬ 
action  potential  Vms{Q)  is  given  in  (2.1)  with  kBT=l.  Notice  that  because  of  the 
orientation  of  imposed  shear  flow,  we  no  longer  have  the  freedom  of  selecting  the 
coordinate  system  to  make  (sin20)  =  0.  Instead,  we  can  select  an  angle  a  such  that 

(sin2(6»-a))  =  0.  (3.2) 

Expanding  (sin26*)  and  (cos2d)  in  (2.1)  yields 

VMs{d)  =  -U{cos2{e-a))cos2{e-a).  (3.3) 

The  potential  caused  by  the  shear  flow  is  a  tilted  periodic  function 

-|0+|l(sR(0)  =  -|0+|-^sin20,  (3.4) 

where  Vsh{&)  =  ^sin2d  is  the  periodic  part  of  the  potential  caused  by  the  shear  flow, 

normalized  by  half  of  the  Peclet  number.  The  case  of  high  shear  rate  (large  Peclet 
number)  was  studied  in  [26,  27]  and  most  recently  analyzed  in  [30].  In  this  work  we 
consider  the  case  where  the  imposed  shear  is  weak  (small  Peclet  number),  which  will 
allow  us  to  do  asymptotic  expansions. 

We  are  interested  in  finding  steady  state  solutions.  We  first  substitute  f{6)  = 
f^^^0)  +  O{e)  into  the  Smoluchowski  equation  (3.1)  and  keep  only  the  leading  order 
terms  to  obtain 

^[<s(^)/^°^  +  ^/^°^]=0,  (3.5) 

where  the  leading  order  term  of  the  Maier-Saupe  interaction  potential  is 

vl^l{e)  =  -U{cos2{9-a))^°'>cos2{e-a).  (3.6) 

Here  the  notation  (cos2(0  — means  the  average  of  cos2(6*  — a)  with  respect  to 
the  leading  order  term  of  the  probability  density  (3.6)  is  simply  the  governing 

equation  for  nematic  polymers  with  only  Maier-Saupe  interaction  potential.  In  par¬ 
ticular,  the  effect  of  shear  does  not  appear  in  the  leading  term  equation  (3.6).  Thus, 
for  0<U  <2,  we  have 

=  vj°U9)  =  o. 


(3.7) 
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For  U  >2,  there  are  two  steady  states  [7,  25]:  the  isotropic  state  and  nematic  state. 
For  U>2,  the  isotropic  state  is  unstable  [11];  the  stable  nematic  state  is 


exp(rcos2(6»-a)) 

J  f27I-  _  , 


/q  ^exp(rcos2(6*  — 


yMU^)  =  -rcos2{9-a), 


(3.8) 


where  r>0  is  the  solution  of  equation  (2.8).  Below  we  are  going  to  continue  with 
the  case  of  17  >  2  (i.e.  the  case  of  r  >  0)  to  see  under  what  condition  a  steady  state 
exists.  Due  to  the  absence  of  the  shear  effect  in  the  equation  for  the  leading  order 
terms,  the  angle  a  in  (3.8)  is  completely  arbitrary  and  cannot  be  determined  from 
the  leading  order  equation  (3.5).  Let  us  try  to  determine  the  angle  a  from  the  full 
equation  instead.  A  solution  of  the  full  Smoluchowski  equation  (3.1)  satisfies 


(3.9) 


where  J{9)  is  the  probability  flux.  At  a  steady  state,  J{9)  =  J  is  a  constant  inde¬ 
pendent  of  9.  Multiplying  both  sides  of  (3.9)  by  the  integrating  factor  exp[—|d-|- 
^VsH{d)  +  VMs{0)]  and  integrating  from  q  to  q  +  2TT,  we  have 


[l-exp(-e7r)]exp(^-|g-h|Vsff(<?)  +  VMs(g))  f{q) 


fq+2Tr 


=  J 


exp 


-2^+  ^Vsh{O)  +  Vms{0) 


d9. 


(3.10) 


Dividing  both  sides  by  J  •exp{—^q+ ^VsH{q)  +  VMs{q)),  expanding  the  right  hand 
side  and  using  the  fact  that  Vms(^*)  =  l^M_g(0) -|-C)(e),  we  obtain 

[1  — exp(— ett) 


J 


-fiq) 


=  exp  ( —V, 


(-4°s(9)) 


exp 


¥^1(9)]  d9-{l  +  0{e)). 


(3.11) 


Integrating  both  sides  from  0  to  27r  and  using  that  f{q)dq=l  gives  us 
[1  — exp(— £7r)] 


J 

exp(-V^°s(g))dg 
=  4Tr'^g{r)-{l  +  0{e)) 

r(0) 


)  •(l  +  0(e)) 


(3.12) 


In  the  above,  we  have  used  =  — rcos2(^  — a)  and  function  g{r)  is  defined  as 


g{r)=^-^  J  exp{rcos29)d9^  . 


(3.13) 


where  r>0  is  the  solution  of  (2.8).  Expanding  1  —  exp(— ett)  =  ett •  (1 -|- 0(e)),  we  ex¬ 
press  the  probability  flux  J  as 


(3.14) 
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To  find  the  condition  for  the  existence  of  steady  state,  we  seek  another  expression  for 
probability  flux  J.  Integrating  both  sides  of  (3.9)  from  0  to  27r,  we  have 

-2  +  2I  VsHiO)f{9)d0  +  J^  VMs{0)f{9)d9  =  -27rJ.  (3.15) 

Let  us  examine  these  terms  one  by  one.  The  third  term  on  the  left  is, 

VMs{9)md9  =  -(-rcos2{9-a))md9 

2rsiii2{9  —  a)f{9)d9  =  2r{sm2{9  —  a))  =  0.  (3.16) 

Here  we  have  used  the  fact  that  in  the  Maier-Saupe  potential,  the  angle  a  is  selected 
to  make  (sin2(6l  — a))  =  0.  This  result  is  not  surprsing.  The  term  V^g{9)f{9)d9 
is  the  total  force  on  the  ensemble  of  polymer  rods  due  to  the  force  from  the  Maier- 
Saupe  interaction  potential.  Since  the  Maier-Saupe  potential  is  caused  by  the  mutual 
interaction  of  polymer  rods  with  each  other,  Newton’s  third  law  guarantees  that  the 
total  force  on  the  ensemble  is  zero.  The  second  term  on  the  left  of  (3.15)  is. 


r27T 


/>27r 


VsH{9)md9=- 


cos29f{9)d9 


£  £ 

=  -(cos26l)  =  -(cos2(0  — a)cos2Q;  — sin2(0  — a)sin2Q;) 


=  ^ (cos 2(0  — a) cos 2a)  =  |  (cos2(0  — a))^°^cos2a-|-0(£) 


-  —cos2a  +  0{e) 


(3.17) 


Here  again  we  used  (sin2(0  — a))  =  0.  Substituting  the  results  (3.16)  and  (3.17)  into 
(3.15)  yields  another  expression  for  the  probability  flux  J 


47r  J  =  e  • 


1—  £cos2a  + 0{s) 


Equating  (3.14)  with  (3.18),  we  obtain 


1 

9{r) 


1 


r 

U 


cos  2a. 


(3.18) 


(3.19) 


(3.19)  is  the  equation  for  determining  the  angle  a.  Here  mathematically  we  treat  r 
as  an  independent  variable  and  treat  others  as  functions  of  r.  In  particular,  C/  is  a 
function  of  r  (although  from  physical  point  of  view,  it  is  more  reasonable  to  treat  U 
as  an  independent  variable  and  treat  r  as  a  function  of  U).  Now  we  investigate  for 
what  values  of  r  equation  (3.19)  has  a  solution  and  for  what  values  of  r  it  does  not. 
Writing  cos2a=  1  —  2sin^a  and  rearranging  terms,  we  arrive  at 

2sin2a=^,  (3.20) 

9[r) 


where  function  h{r)  is  defined  as 


h{r)  = 


U  r 


l-6f(r)(l 


r 


(3.21) 


H.  ZHOU  AND  H.  WANG 


119 


The  properties  of  h{r)  and  g{r)  are  described  in  the  theorem  below. 

Theorem  3.1.  For  function  g{r)  defined  in  (3.13)  and  function  h{r)  defined  in 
(3.21),  we  have 

1.  (jf(0)  =  l  and  g{r)  >0. 

2.  a  (r)  >  0  for  r  >  0,  whieh  implies  that  q(r)  >  1  for  r  >  0. 

3.  /i(0)  =  l  and  /i'(0)  =  -l<0. 

4.  lim  h(r)  =  —00. 

r — >-(x> 

The  proof  of  Theorem  3.1  is  given  in  Appendix  A. 


Fig.  3.1.  Graphs  of  functions  h{r)  and  The  left  panel  shows  the  graph  of  function  h{r). 

From  the  graph,  it  is  clear  that  h{r)  is  a  decreasing  function  of  r  for  r>0  and  h{r)—^—oo  as 
r— »'+oo,  as  described  in  Theorem  1.  h{r)  intersects  with  the  horizontal  axis  at  ro  =  1.32757384. 
Function  h(r)  satisfies  l>h(r)>0  for  r<ro  and  h{r)<0  for  r>rQ.  The  right  panel  shows  the 
graph  of  function  . 


The  graph  on  the  left  in  Figure  3.1  indicates  that  h{r)  is  a  decreasing  function  of 
r  for  r>0.  A  rigorous  mathematical  proof  for  this  conclusion,  however,  is  still  open. 
Since  h(0)  =  l,  h{r)  is  a  decreasing  function  of  r  for  r>0  and  lim  h(r)  =  — oo,  there 

r— ^oo 

exists  a  unique  positive  number  tq  such  that 


1  >  h{r)  >  0 

for 

r  <ro, 

h{r)  <0 

for 

r  >ro. 

Using  the  property  that  (7(0)  =  !  and  g{r)>l  for  r>0,  we  see  that  the  right  side  of 
equation  (3.20)  has  the  behavior 

h(r) 

l>4^>0 

9h) 

for 

r  <ro, 

h(r) 

)  :<o 

for 

r  >ro. 

This  behavior  of  is  shown  in  the  graph  on  the  right  in  Figure  3.1.  Let  Uq  =  U{ro) 
be  the  value  of  U  corresponding  to  r  =  ro  where  function  U{r)  is  defined  in  (2.10). 
Since  U{r)  is  an  increasing  function  of  r,  the  region  of  U  <Uq  corresponds  to  the 
region  of  r  Kvq. 
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Now  let  us  get  back  to  solving  equation  (3.20).  From  the  behavior  of  obtained 
above,  it  is  clear  that  when  U  <  Uq,  a  steady  state  solution  exists.  When  U  >  Uq,  there 
is  no  steady  state  solution.  The  numerical  values  of  tq,  Uq  and  the  corresponding  order 
parameter  can  be  calculated  using  numerical  integration  and  are  listed  below: 


ro  =  1.32757384, 

C/o  =  C/(ro)  =  2.41144646, 

so  =  ^  =0.550530090.  (3.22) 

Uo 

Our  value  of  Uq  is  consistent  with  the  one  found  in  [26,  27]. 

4.  Multi-scale  solution  of  the  Smoluchowski  equation  for  nematic  poly¬ 
mers  under  a  weak  shear 

The  analysis  in  the  previous  section  shows  that  under  an  imposed  weak  shear 
when  the  normalized  polymer  concentration  U  is  above  the  threshold  Uq  there  is 
no  steady  state  solution.  To  study  the  unsteady  state  dynamics  of  the  probability 
density  for  the  case  of  U  >Uo  and  to  study  the  convergence  of  the  probability  density 
to  steady  state  for  the  case  of  U  <Uo,  we  carry  out  multi-scale  asymptotic  expansion 
[2,  21]  of  the  Smoluchowski  equation  (3.1).  The  approach  of  multi-scale  expansions 
has  been  applied  to  study  the  nonlinear  dynamics  of  a  nematic  liquid  crystal  under  a 
shear  flow  for  a  Landau-de  Gennes  model  [28]  and  a  Doi-Hess  mesoscopic  orientation 
tensor  model  [3] .  Our  goal  here  is  to  capture  the  slow  time  evolution  of  the  probability 
density  governed  by  the  Smoluchowski  equation. 

Let  To  =  t  and  Ti=e-t.  We  seek  an  expansion  of  the  form 

/(0,t)  =  /(°)(0,To,Ti)  +  £/W(0,ro,ri)  +  O(e2).  (4.1) 


Substituting  the  multi-scale  expansion  (4.1)  into  the  Smoluchowski  equation  (3.1)  and 
keeping  only  the  0(1)  terms,  we  have 


df(o) 


_9 


w(o)'  f(o) 

^MS  J 


IL  fio) 

de^ 


(4.2) 


This  is  exactly  the  governing  equation  for  nematic  polymers  in  the  absence  of  shear 
flow.  Steady  state  solutions  of  (4.2)  are  equilibria  and  are  given  by  the  Boltzmann 
distribution 

/p  exp (r COS 2 (6*  — a)) 

Again,  mathematically  it  is  more  convenient  if  we  simply  treat  r  as  an  independent 
variable.  If  we  choose  to  treat  U  as  the  independent  variable  and  view  r  as  a  function 
of  U,  then  r{U)  can  be  solved  from  equation  (2.8).  Suppose  the  polymer  orientation 
starts  with  the  equilibrium  distribution  given  in  (4.3).  Equation  (4.2)  dictates  that 
the  polymer  orientation  distribution  stays  unchanged  with  respect  to  Tq.  That  is. 


a/(o)(0,To,ri) 

dTo 


(4.4) 


We  say  that  f^^\9,To,Ti)  is  a  quasi-steady  state  if  it  is  independent  of  Tq.  In  particu¬ 
lar,  the  phase  angle  a  in  (4.3)  is  independent  of  Tq.  However,  as  we  will  see,  the  phase 
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angle  a  may  evolve  with  Ti,  which  corresponds  to  the  slow  time  evolution  caused  by 
the  weak  shear  flow.  To  calculate  a{Ti),  we  look  at  the  expansion  of 

df{e,t)  _df^°H0,To,T,)  p/(o)(g,To,ri)  9/(i)(g,ro,ri)i 

dt  dTo  dTi  dTo  \ 

If  the  polymers  start  with  an  equilibrium  distribution  for  pure  nematic  polymers  (no 
shear),  then  /(‘’^(0,To,Ti)  is  a  quasi-steady  state  and  the  0(1)  term  in  the  expansion  of 
is  zero.  The  0(e)  term  in  contains  two  parts.  In  multi-scale  expansions, 

to  separate  the  effects  of  time  scales  Tq  and  Ti,  it  is  necessary  to  impose  condition(s) 
on  the  /(^^(0,To,Ti)  term.  Here  we  impose  the  condition  that  f^^\9,TQ,Ti)  is  also  a 
quasi-steady  state  (that  is,  it  is  independent  of  Tq).  From  physical  point  of  view,  it  is 
reasonable  to  impose  the  condition  that  f^^^6,To,Ti)  is  a  quasi-steady  state.  In  the 
absence  of  shear  flow,  the  Maier-Saupe  interaction  potential  drives  the  polymers  to 
the  equilibrium  orientation  distribution.  The  effect  of  imposed  weak  shear  flow  will 
disturb  the  equilibrium  distribution.  Since  the  imposed  shear  is  weak  the  perturbation 
on  orientation  distribution  is  also  small.  More  specifically,  during  a  time  period  of 
T,  the  total  perturbation  on  the  orientation  distribution  is  of  the  order  er.  During 
the  time  period  of  r,  the  orientation  distribution  will  constantly  try  to  relax  to  new 
equilibriums.  As  we  discussed  above,  the  leading  term  of  the  orientation  distribution 
is  in  a  quasi-steady  state  (independent  of  Tq).  The  leading  term  of  the  orientation 
distribution  rotates  with  respect  to  the  slow  time  scale  Ti.  In  a  time  period  of  t, 
the  leading  term  of  the  orientation  distribution  rotates  by  £r.  Here  it  is  important  to 
point  out  that  the  Maier-Saupe  interaction  potential  rotates  along  with  the  orientation 
distribution.  So  if  we  rotate  the  coordinate  system  to  follow  the  leading  term  of 
the  orientation  distribution,  then  the  leading  term  of  the  Maier-Saupe  potential  will 
be  stationary.  In  the  rotating  coordinate  system,  the  effect  of  the  weak  shear  also 
changes  at  most  by  0(£- er)  =  0(£^r).  As  a  result,  in  the  rotating  coordinate  system, 
the  orientation  distribution  can  change  at  most  by  0(£^r)  over  the  time  period  of 
T.  Therefore,  in  the  rotating  coordinate  system,  the  rate  of  change  of  the  orientation 
distribution  is  of  the  order  0{e^),  which  implies  that  f^^\6,To,Ti)  is  independent 
of  the  fast  time  Tq.  It  has  also  been  verified  in  our  extensive  numerical  simulations 
(results  not  shown)  that  /^^^(0,To,Ti)  is  a  quasi-steady  state. 

Since  f^^'>{9,TQ,Ti)  is  independent  of  Tq,  the  expansion  of  becomes 

df{9,t)  _a/(o)(0,To,Ti)  2^ 

dt  dTi  ^  ^ 

Using  the  Boltzmann  form  of  /*^°^(6*,To,Ti)  given  in  (4.3),  we  have 

df{9,t)  Jf(^H9,TQ,T,)  dajn)  , 

dt  d0  dTi  ''  ^  ^ 

Substituting  the  expansion  of  into  the  Smoluchowski  equation  (3.1),  moving 

all  terms  to  the  right  side  and  using  the  fact  that  a  is  independent  of  9  (a  depends 
only  on  Ti),  we  obtain 

O=^l^(-l  +  lv'^(O)  +  V;,smf(9)  +  ^+ef«^H0)^^+O(s^)'j.  (4.6) 

It  follows  that 

(- 1  +  lv'M9)  +  Ums(0))  m+^+  (^) +  O(s^)  =  -eh  (4.7) 
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where  e  J2  is  the  0{e)  term  of  the  integration  constant,  independent  of  6.  Since  the 
solution  of  the  0(1)  term  equation  (4.2)  is  an  equilibrium,  the  0(1)  term  of  the 
integration  constant  is  zero. 

In  (4.7),  both  and  J2  are  unknown.  We  need  to  solve  for  both  of  them 

simultaneously.  For  that  purpose,  we  need  two  linear  equations  for  and  J2. 

Integrating  (4.7)  from  0  to  27r  and  dividing  by  e  yields 


-  i  i  ^  cos(2a(ri ) )  +  +  0(e)  =  -2nJ^ . 


(4.8) 


In  the  integration  above  we  have  used  results  (3.16)  and  (3.17).  The  second  linear 
equation  for  and  J2  is  obtained  by  dividing  (4.7)  by  ef{6),  approximating  f{6) 

by  and  then  integrating  from  0  and  27r: 

da{Ti) 


-tv  +  2tt  +0{e)  =  -(2Trfg{r)J2. 

dl  1 


(4.9) 


Here  we  have  used  the  fact  that 


/(O)(0) 


d0  = 


exp{rcos2{9  —  a))d9 


exp{—rcos2{9  —  a))d9 


=  (27r)^5(r).  (4.10) 

where  function  g{r)  is  defined  in  (3.13).  Solving  from  (4.8)  and  (4.9)  yields 


da{Ti) 

dTi 


1- 


9{r) 

g{r)-l  U 


-cos(2a(ri)) 


Writing  cos(2a(ri))  =  1  —  2sin^Q;(Ti),  we  obtain 


da{Ti) 

dTi 


gjr)  r 
g{r)-l  U 


sin2(a(Ti))- 


(4.11) 


(4.12) 


where  function  h{r)  is  defined  in  (3.21).  The  slow  time  evolution  of  the  phase  angle 
Q;(ri)  is  governed  by  the  differential  equation  (4.12).  In  equation  (4.12),  the  factor 
■  xj  is  always  positive  (for  U  >2).  As  we  discussed  in  the  previous  section,  for 
U  <Uo  =  2.41144646,  we  have  0  <  <  0.5  and  q;(Ti)  will  converge  to  a  steady  state 

angle  satisfying 


Mr) ' 


(4.13) 


There  are  two  angles  satisfying  (4.13)  in  the  interval  [— f,f],  one  positive  and  one 
negative.  The  whole  problem  is  periodic  with  period  tt  so  we  only  need  to  examine 
an  interval  of  size  tt.  The  stability  of  a  steady  state  angle  is  determined  by 
expanding  equation  (4.12)  around  Near  we  have 


d{a{Ti)  — 

Jfi 


9{r) 

9{r)-l 


■  ^  •sin(2Q;*''®^)(a(Ti)  —  Qf^'®'^). 


It  follows  that  is  stable  if  sin(2a(^^)  <0  and  is  unstable  if  sin(2a(‘^))  >0.  So 
the  positive  steady  state  angle  is  unstable  whereas  the  negative  steady  state  angle  is 
stable. 
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When  U>Uo  =  2.41144646,  we  have  <  0  and  q;(Ti)  will  keep  increasing  with 
Ti  because  the  right  side  of  (4.12)  is  always  positive.  So  when  U  >  t/g  =  2.41144646, 
the  director  will  keep  rotating  (tumbling)  and  there  is  no  steady  state.  In  the  following 
we  derive  the  period  of  the  tumbling  state. 

Let  Tp  be  the  period  of  tumbling  in  slow  time  Ti  =  et.  For  U  >  Uo  =  2.41144646, 
differential  equation  (4.12)  has  the  form 


da 


Cl  (sin^a  +  C2) ,  ci  = 


_  9{r)  r  _-Hr) 

•  77  >  U,  C2  =  „  ,  ,  >  0. 


g{r)-l  U 


(4.14) 


Consider  the  period  [tq,tq  +  Tp\  in  Ti  that  corresponds  to  a  going  from  0  to  27r. 
Dividing  both  sides  of  (4.14)  by  (sin^a  +  C2)  and  integrating  over  the  period  [ro,To  +  rp] 
yields 


/.27I- 

/  — 2 - da  =  ciTp.  (4-15) 

Jo  sin  a  +  C2 

Because  of  the  symmetry  of  function  sin^a,  on  the  left  side  of  (4.15),  the  integral  over 
[0,27r]  is  equal  to  4  times  the  integral  over  [0,7r/2].  Using  substitution  s  =  tan(a),  we 
obtain 


1 


/o  (1  +  C2)s2  +  C2 


ds  =  CiTp. 


(4.16) 


Applying  substitution  u  =  s  and  using  the  fact  that  ^2^1  du  =  f,  we 

arrive  at 

rp  =  27T—^^==.  (4.17) 

Cl  VC2(1  +  C2) 


Fig.  4.1.  Solid  line  shows  the  normalized  tumbling  frequency  (kje)  as  a  function  of  U.  The 
dashed  line  shows  the  asymptotic  approximation  for  the  normalized  tumbling  frequency  as  U  — >oo. 


Thus,  in  time  the  period  of  tumbling  is  given  by 


t 


p  — 


27r 

s 


(4.18) 
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The  tumbling  frequency  (fc)  normalized  by  the  Peclet  number  (e)  is 

k  I  I  f  g{r)  r\  j  / -h{r)\  f 

£  tpE  27r  \g{r)-l  u)\l\2g{r)j[  \2g{r)  J 


(4.19) 


In  Figure  4.1,  the  solid  line  shows  the  normalized  tumbling  frequency  (fc/e)  as  a 
function  of  U.  As  U  increases  from  Uq,  the  normalized  tumbling  frequency  increases 
steeply  to  reach  a  maximum  value.  Beyond  the  location  where  the  maximum  value 
is  attained,  when  U  increases  further  to  infinity,  the  normalized  tumbling  frequency 
decreases  monotonically  to  zero.  As  U^oo,  the  order  parameter  satisfies  (cos  20)  = 
r/U—i'l,  which  implies  r—>-oo.  In  Appendix  A,  we  show  that  as  r^oo  functions 

g(r)  and  h{r)  have  asymptotic  expansions  g(r)  =  l/(27rr)exp(2r)(lH - )  and  h{r)  = 

—  l/(47rr^)exp(2r)(H - ).  Substituting  these  results  into  (4.19),  we  arrive  at  an 

asymptotic  approximation  for  the  normalized  tumbling  frequency  as  U  —>■  oo. 

fell 

-  = - =H -  ast/^+oo.  (4.20) 

e  4tt 


In  Figure  4.1,  the  dashed  line  shows  the  asymptotic  approximation  for  the  normalized 
tumbling  frequency  as  U  oo.  As  shown  in  the  figure,  the  asymptotic  approximation 
is  very  close  to  the  true  value  for  U  >4. 

In  summary,  the  leading  term  expansion  of  the  solution  of  Smoluchowski  equation 
(3.1)  is  given  by 

/(°)(0.t)=  yxp|.cos2(0-»(rt))]  _ 

e:xjp[r  cos2{0  —  a{et))]d9 

where  the  anglular  coordinate  of  the  director  q;(Ti)  is  governed  by  the  differential 
equation  (4.12).  For  U  <Uo,  the  director  converges  to  a  stationary  angle  and  we 
have  a  steady  state.  For  U  >  Uq,  the  director  keeps  rotating  and  we  have  a  tumbling 
state.  For  2D  nematic  polymers  under  an  imposed  weak  shear,  these  are  the  only  two 
possible  states.  In  particular,  there  is  no  transient  oscillatory  (wagging)  solution. 

5.  Behavior  of  the  Smoluchowski  equation  near  U  =  2 

As  discussed  in  Section  2,  in  the  absence  of  shear  flow,  the  isotropic  to  nematic 
phase  transition  occurs  at  U  =  2.  For  fixed  U  <2,  as  e^O,  the  polymer  orientation 
distribution  converges  to  the  isotropic  phase.  For  fixed  C/>  2,  as  0,  the  orientation 
distribution  converges  to  the  nematic  phase.  Now  we  study  the  behavior  near  U  =  2 
when  £  is  small  but  is  finite. 

Notice,  in  particular,  that  in  the  region  near  U  =  2  we  have  U  <Uo  and  conse¬ 
quently  the  polymer  orientation  distribution  has  a  steady  state.  Here  again,  it  is 
mathematically  more  convenient  to  treat  r  as  the  independent  variable.  Below  we 
investigate  the  relation  between  U  and  r  near  r  =  0  in  the  steady  state.  More  specifi¬ 
cally,  we  study  the  region  where  r  is  the  range  of  [0,  o(£s)].  In  the  proof  of  Theorem 
1  given  in  Appendix  A,  we  showed  ^(0)  = /i(0)  =  1,  which  implies  limr^o  2g(r) 

It  follows  that  when  r  is  small,  the  stable  steady  state  phase  angle  determined  by 
equation  (4.13)  is  approximately  a  =  —\  (the  unstable  steady  state  phase  angle  is  ap¬ 
proximately  a=  ^).  Substituting  the  phase  angle  a  into  the  Maier-Saupe  interaction 
potential,  we  have 

TT 

—r  cos  2(0  — a)  =  — rcos(20-|-  — )  =  rsin20. 
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That  is,  when  r  is  small,  the  periodic  part  of  the  potential  caused  by  the  shear  flow 
is  aligned  with  the  Maier-Saupe  potential.  The  sum  of  these  two  potentials  is 

^-VsH{e)  +  VMs{0)  =  -{^-+r)cos2{9-a).  (5.1) 

If  we  treat  r  as  an  independent  variable,  [/  as  a  function  of  r  is  given  by 

{cos2{e-a)); 

The  subscript  in  [4(r)  indicates  that  it  is  for  the  case  where  the  shear  flow  is  present 
and  the  effect  of  shear  flow  is  included  in  calculating  the  average.  The  inverse  function 
rs{U)  is  probably  more  meaningful  from  the  physical  point  of  view  and  the  order 
parameter  s  as  a  function  of  the  normalized  polymer  concentration  U  is  given  by 
s{U)  =  .  For  small  r,  it  is  necessary  to  include  the  effect  of  shear  flow  in  the 

leading  order  expansion.  In  (5.1),  when  r  is  small,  e  may  not  be  ignored  and  may 
even  be  the  dominant  term  relative  to  r.  To  include  the  effect  of  shear  flow,  we  start 
by  finding  a  more  accurate  expansion  for  the  steady  state  probability  density  f{0). 
Note  that  in  the  presence  of  shear  flow,  the  steady  state  probability  density  is  no 
longer  given  by  the  Boltzmann  distribution.  Fortunately,  f{6)  can  be  expressed  from 
equation  (3.10)  in  the  form: 

lo^ exp  [- fu;  +  I VsH {oj  +  9)  +  Vms{‘^  +  0)]  duj 
e^p{lVsH{O)  +  VMs{0)) 


Here  we  have  absorbed  into  the  constant  c  all  factors  that  are  independent  of  9.  Notice 
that  (3.10)  is  exact  so  (5.2)  is  also  exact.  Now  let  us  expand  the  integral  in  (5.2)  and 
keep  terms  up  to  0{e).  Remember  that  r  is  in  the  range  of  [0,  o(eJ)]. 


2tt 


exp 

2-k 

exp 

2tt  , 


—  -LO+  -V5/i-(a;  +  6*)  +  Vms(‘^  +  ^) 


doj 


£  £ 

— -w  —  (-  +  r)  cos2{uj  +  9  —  a) 


duj 


£  £  £ 

1  —  -uj  —  e{-  +  r)  u}COs2{uj  +  9  —  a)  —  (  -  +  r)^cos^2(w  +  0 


cr)  + 


duj. 


In  the  above,  only  the  third  term  may  depend  on  9  but  it  is  of  the  order  e-0(|  + 
r)  =  o{ei).  All  other  terms  are  independent  of  9.  In  particular,  the  fourth  term  is 
independent  of  9  and  is  of  the  order  0((| +  r)^)  =  o(£5 ).  Thus,  accurate  up  to  o{e^) 
the  integral  in  (5.2)  is  independent  of  9.  In  (5.2),  absorbing  the  integral  into  the 
constant  c,  we  obtain 


f{9)  =  c-exp(^-^VsHi0)  -  Vms{0))  +o(e4) 

_  exp[(|+r)cos2(6>-a)]  5 

fg^exp  (|  +  r)cos2(0  — a)  d9 


Comparing  the  probability  density  in  (5.3)  with  the  probability  density  in  the  case 
of  no  shear  flow  given  in  (2.6),  we  see  that  the  only  difference  is  that  r  in  (2.6)  is 
replaced  by  (|  +  r)  in  (5.3).  As  a  result,  the  order  parameter  in  the  presence  of  shear 
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flow  is  equal  to  the  order  parameter  in  the  absence  of  shear  flow  but  with  r  replaced 
by  ( I  +  r)  when  calculating  the  average.  Mathematically,  we  have 

{cos2{0-a))s=  (cos2(6»-a))|^^(j_^^).  (5.4) 

Using  =  {cos2{9  —  a)) e  and  =  (cos2(0  — a)),  it  follows  that 

r  _  l+r 
Ue{r)  C/(f  +  r)' 

Writing  Us{r)  as  a  function  of  r,  we  arrive  at 

t/,(r)  =  ^.C/(|+r) 

where  function  Ug{r)  is  for  the  case  with  shear  flow  and  function  U{r)  is  for  the  case 
of  no  shear  flow.  Relation  (5.6)  shows  that  the  function  Ueir)  can  be  obtained  by  a 
simple  algebraic  transform  on  the  function  U{r).  Relation  (5.6)  is  also  valid  when  r  is 
not  small.  When  r  is  finite  and  as  0,  (5.6)  simply  tells  us  that  t4(r)  ~U{r),  which 
is  true  and  is  what  we  obtained  in  our  leading  term  expansion.  Therefore,  relation 
(5.6)  is  valid  for  the  full  range  of  r. 

The  phase  diagram  can  be  expressed  in  terms  of  the  inverse  function  of  Ug{r): 


(5.5) 

(5.6) 


Ss{U) 


_  rejU) 
U 


Phase  diagrams  for  £  =  0  (no  shear)  and  e  =  0.1  (weak  shear)  are  shown  in  Figure  5.1. 
There  are  several  regimes  in  the  phase  diagram.  To  analyze  these  regimes,  we  need 
the  expansion  of  function  U{r)  around  r  =  0: 


U{r)  =  2 


(^l+lr^  +  0{r^)y 


(5.7) 


Regime  1:  r  is  in  the  range  of  [0,  O(-ye)]. 

In  this  case,  expanding  {7(|  +  r)  in  (5.6)  yields 


Ue{r)  = 


■2{l  +  0{e)), 


(5.8) 


which  indicates  that  the  corresponding  range  for  Ug  is  [0,  2  —  0{y/e)].  In  (5.8),  ne¬ 
glecting  the  0{s)  term  and  expressing  r  as  a  function  of  [4,  we  have 


U 


4  2-U’ 


TejU) 

u 


4  2-U 


(5.9) 


(5.9)  is  valid  for  U  in  the  range  of  [0,  2  —  0{y^)].  In  particular,  we  have  Se(0)  =  |. 
In  other  words,  the  order  parameter  is  non-zero  even  when  the  normalized  polymer 
concentration  U  goes  to  zero.  Here  the  alignment  is  caused  by  the  imposed  shear 
flow. 

Regime  2:  r  is  of  the  order  of  0(e^/^). 

In  this  case,  expanding  U{^  +  r)  in  (5.6)  yields 


Ue{r)  =  - 


1  +  xr^- 


4 


(5.10) 


H.  ZHOU  AND  H.  WANG 


127 


U  U 


Fig.  5.1.  Phase  diagrams  for  the  2D  nematic  polymers  with  and  without  shear  flow,  e  is  the 
Peclet  number.  The  left  panel  shows  r  as  a  function  of  the  normalized  polymer  concentration  U 
where  r  is  defined  as  r  =  U{cos29).  The  right  panel  shows  the  order  parameter  s  as  a  function  of 
U  where  the  order  parameter  is  defined  as  s  =  {cos20).  In  the  absence  of  shear  flow,  the  isotropic- 
nematic  phase  transition  occurs  atU  =  2.  In  the  presence  of  shear  flow,  the  phase  transition  occurs 
gradually  over  a  region  near  U  =  2. 


which  indicates  that  the  corresponding  range  for  C4  is  [2  — 2  +  0(£^/^)]  and 
the  corresponding  order  parameter  is  of  the  order  of  From  (5.10),  the  order 

parameter  s  can  be  related  to  the  normalized  polymer  concentration  U  through  a 
cubic  equation.  In  (5.10),  dividing  both  sides  by  U  and  using  r  =  U  -s  yields 


1  = 


l  +  U-s 


Using  s  =  0(e^/^)  and  17  =  2  +  0(e^/^)  to  write  1+  =  1+  +  0(£‘^/^),  neglecting 

the  term,  and  then  solving  U  in  terms  of  s,  we  arrive  at 


t/s(s)  =  2+( 


(5,11) 


The  inverse  function  s^{U)  can  be  obtained  by  solving  cubic  equation  (5.11): 


Se 


where  the  quantity  q  is 


(5.12) 


/ley/^  u-2 


The  value  of  s  that  corresponds  to  C/  =  2  is 


128 


2-D  NEMATIC  POLYMERS  UNDER  AN  IMPOSED  SHEAR 


As  we  pointed  out  earlier  in  the  discussion,  (5.6)  is  the  relation  for  the  stable 
solution.  For  the  unstable  solution,  we  have  a  similar  relation 

U,{r)  =  ^-U{r-j),  (5.14) 

r  4  4 

where  Ue{r)  corresponds  to  the  unstable  solution  in  the  presence  of  shear  flow  and 
function  U (r)  is  for  the  nematic  solution  with  no  shear  flow.  Below  we  show  that  the 
unstable  branch  Ug{r)  has  a  minimum.  When  r  is  of  the  order  of  0(e^/^),  expanding 
(5.14)  and  neglecting  0(£‘*/^)  terms  yields 


U,{r)-2=^  +  ^r^.  (5.15) 

In  this  case,  the  two  terms  on  the  right  hand  side  of  (5.15)  are  of  the  same  order  and 
[4  is  of  the  order  of  2  +  0(£^/^).  In  (5.15),  minimizing  [4  with  respect  to  r,  we  obtain 
that  C4  attains  the  minimum  at  r  =  £^/^.  The  minimum  value  of  [4  (th®  unstable 
branch)  is 


min[/e  =  2+^£^/^.  (5.16) 

For  r«  the  first  term  on  the  right-hand  side  of  (5.15)  is  the  dominant  term 

while  for  the  second  term  is  the  dominant  term.  For  both  of  these 

two  caess,  we  have  Ug  —  2>>0{e^^^).  Therefore,  the  minimum  of  C4  attained  at 
r  =  £^/^  is  the  global  minimum  of  C4-  The  scaling  relation  (5.16)  is  consistent  with 
the  one  obtained  in  [23]  for  the  Doi  closure.  Phase  diagrams  of  the  stable  branch, 
the  unstable  branch  and  the  travelling  wave  branch  (unsteady,  periodic  in  time)  for 
£  =  0.1  are  shown  in  Figure  5.2. 

6.  Conclusion 

In  this  paper  we  have  studied  analytically  the  2D  Smoluchowski  equation  under 
an  imposed  weak  shear.  Our  main  results  can  be  summarized  as  follows. 

•  In  the  presence  of  an  imposed  weak  shear  there  exists  a  critical  value 
([/q  =  2.41144646)  for  the  normalized  polymer  concentration  U.  Steady  state 
solution  exists  only  for  U  <Uq. 

•  The  multi-scale  asymptotic  analysis  describes  the  slow  time  evolution  driven 
by  the  weak  shear.  It  is  found  that  the  order  parameter  of  the  orientational 
distribution  of  rodlike  molecules  is  time-invariant  while  the  director  changes 
with  time.  It  is  further  shown  that  when  U  <Uo,  the  director  of  the  orienta¬ 
tional  distribution  converges  to  a  stable  steady  state  position;  when  U>Uo, 
the  angular  velocity  of  the  director  is  position-dependent  and  is  always  pos¬ 
itive.  As  a  result,  the  orientational  distribution  is  temporally  periodic  and 
can  never  reach  a  steady  state  for  U>Uq. 

•  Finally,  it  is  revealed  that  the  phase  relation  under  weak  shear  can  be  obtained 
from  the  phase  relation  with  no  shear  by  a  simple  algebraic  transformation. 


Appendix  A.  Proof  of  Theorem  3.1. 

Below  we  prove  the  4  items  in  Theorem  1  one  by  one. 

1)  It  is  straightforward  to  verify  that  5(0)  =  !  and  g{r)  >0  from  (3.13). 
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Fig.  5.2.  Phase  diagrams  of  the  stable  steady  state,  unstable  steady  state  and  the  traveling  wave 
(unsteady)  for  the  2D  nematic  polymers  with  £  =  0.1.  Here  e  is  the  Peclet  number.  The  left  panel 
shows  r  as  a  function  of  the  normalized  polymer  concentration  U  where  r  is  defined  as  r  =  U{cos26). 
The  right  panel  shows  the  order  parameter  s  as  a  function  of  U  where  the  order  parameter  is  defined 
as  s  =  (cos26).  Unstable  steady  state  does  not  exist  for  U  <2+|£^/^.  For  U>2-\-^e^^^,  unstable 
steady  states  appear.  A  unstable  steady  state,  when  perturbed,  will  either  converge  to  a  stable  steady 
state  or  converge  to  a  travelling  wave  (unsteady,  periodic  in  time).  The  lower  half  of  the  unstable 
steady  branch  extends  to  U  =  +oo,  which  corresponds  to  the  unstable  isotropic  state  in  the  absence  of 
shear  flow.  The  upper  half  of  the  un  stable  steady  branch  stops  at  Uq  =  2.41144646.  The  stable  steady 
state  branch  also  stops  at  Uq  =  2.41144646.  For  U  >Uo  =  2.41144646,  there  is  no  steady  nematic  state 
(stable  or  unstable).  Instead,  there  is  a  travelling  wave  nematic  state  (unsteady,  periodic  in  time). 


2)  Differentiating  (3.13),  we  have 


dr 


g{r)  =  2^-T  J  exp{rcos26)d9 
=  2^—  J  exp{rcos29)d6 


—  2g{r){cos29). 


J  cos29exp{rcos29))d9 

^  fficos29exp(rcos29)d9 
ffi  exp(rcos29))d9 


(A.l) 


Thus,  to  show  that  g  (r)  >  0  for  r  >  0,  we  only  need  to  show  (cos 20)  >  0  for  r  >  0.  Using 
the  probability  density  given  in  (2.6),  it  is  straightforward  to  verify  that  (cos 20) ^=0  = 
0.  Differentiating  (cos  20)  with  respect  to  r,  we  get 


A(cos20)  =  (cos20(cos20  —  (cos20)))  =  var(cos20)  >  0,  (A. 2) 

dr 

where  var(cos20)  denotes  the  variance  of  random  variable  cos 20.  It  follows  that 
(cos20)  is  an  increasing  function.  In  particular,  we  have  (cos20)  >0  for  r  >0. 

3)  Using  equation  (2.7),  we  write  function  h(r)  in  the  form 

^  (cos20)  [^~g(O(l-(cos20))].  (A.3) 


Using  the  probability  density  given  in  (2.6),  it  is  straightforward  to  verify  that 
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(cos20)r=o  =  O,  {cos^20)r=o=  \  and  (cos^20)r=o  =  O.  Furthermore,  we  have 


(cos  20) 
dr 


r—0 


=  ((cos^  29)  —  (cos26*)^)  “  2  ^ 


=  (cos^20(cos20—  (cos20)))|r=o 


r—0 


— 2(cos20)|r=o  ■  (cos20(cos20—  (cos20)))|r=o  =  0; 
=  2g(r)(cos20)|r=o  =  O, 


r=0 


r—0 


2-^g{r){cos29)  +  2g{r)^{cos26) 


=  1. 


r=0 


(A.4) 


(A.5) 


Using  these  results  to  expand  (cos20)  and  g{r)  around  r  =  0,  we  have 

(cos20)  =  ^r  +  0(r^) 

0(L)  =  l+y +  0(r3). 

Substituting  these  expansions  into  (A. 3),  we  obtain 

h{r)  =  1  — r  +  0(r^). 

Thus,  we  have  h  (0)  =  — 1  <  0. 

4)  In  order  to  find  out  the  asymptotic  behavior  of  h{r)  as  r^+oo,  we  use  Watson’s 
lemma  [2].  First,  we  make  a  change  of  variable: 

<7  =  — cos20+l.  (A. 6) 

As  r^+oo,  the  dominant  contribution  from  integral  exp{rcos29)dd  comes  from 
the  neighborhood  near  cos  20=1  (that  is,  the  neighborhood  near  <7  =  0).  Differentiat¬ 
ing  (A. 6)  and  expanding  for  small  q,  we  have 


2d9  = 


Substituting  into  integral  exp (r  cos  29)  d9  yields 

exp(rcos29)d9  =  4  exp(rcos29)d9 

/■^  11 

=  2exp(r) J  exp{-rq)-^{l  +  -q+---)dq 

/•2  ill 

=  V2exp{r)  exp{-rq){q~^  + -q^ -\ - )dq 

Jo  4 


^=d9=^0  +  lq+---)dq. 


1- 


=  V2exp{r) 


r(^)  ,  ir(i) 
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Similarly,  we  have 


cos  20  exp(r  cos  29)d9  =  A  cos  20  exp(r  cos  29)d9 

/■^  11 
=  2exp(r) J  exp{-rq){l-q)-^{l+-q+---)dq 

r‘2.  ^  3  1 

=  \/2exp(r)  /  exp(-r(7)(g“5_  g,2-| - '^dq 

Jo  4 


=  C2»pMi''<h  3r(i) 


7*  2 


/  W-,  3 

=  ^exp(r)(l-— +  •••) 
Vr  8r 

The  above  two  expansions  lead  to 


as  r^+oo. 


(cos  20)  = 


cos  29exp  (r  cos  29)  d9 

fg^  exp  (r  cos  29)  d9 

1—  — H 1 

=l-±  +  ... 
l+i  +  -  2r 


as  r^+oo 


and 


g(r)=^-^  J  exp(rcos20)(i0^  = 

=  A— exp(2r)(l  +  ^  + 

27rr 


1  2 


^=exp(r)(l  +  — +  •••) 
V27rr  or 

as  r^+oo. 


Substituting  the  expansions  for  (cos 20)  and  g(r)  into  h{r),  we  arrive  at 

1 


h{r)  = 


(cos  20) 
1 


[l-5W(l-(cos20))] 


exp(2r)  - 


as  r^+oo. 


This  completes  the  proof  of  Theorem  3.1. 


(A.7) 


(A.8) 


(A.9) 
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